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ABSTRACT 


To  assist  propeller  designers  in  developing  blade  designs  based  on  composite 
structures,  an  overview  of  composite  material  characteristics  and  available  structural 

•  analysis  methods  tvas  undertaken.  Using  detailed  analyses  of  a  thick-shell  composite  and 
a  solid  bronze  blade,  it  was  shown  that  three-dimensional  finite-element  analysis  is  useful 
for  predicting  relevant  stresses  and  deflections.  A  more  simplified  approach  was  not 

•  satisfactory.  To  establish  the  feasibility  composite  blades,  further  information  about 
acceptable  stress  and  deflection  levels  is  needed. 


ADMINISTRATIVE  INFORMATION 

This  report  is  submitted  in  paitial  fulfillment  of  Task  2  of  Project  RB23C22, 
Program  ND3A/PE0602323N.  The  wort  described  herein  was  sponsored  by  the 
Office  of  Naval  Technology  (ONT  233)  and  performed  by  the  David  Taylor 
Research  Center  (DTRC)  under  Work  Units  1506-060  and  1506-160,  DN501143. 

INTRODUCTION 

Composite  materials  are  defined  as  the  man-made  products 
manufacture  by  the  assembly  of  at  least  two  dffierent  materials  that  are 
chemically  distinct  on  a  macroscopic  scale  and  have  a  clearly  recognizable 
interfiice  between  them.l^*21  A  typical  construction  involves  a  reinforcing  fiber 
agent  embedded  in  a  compatible  matrix  material  such  as  a  resin.  In  general, 
composite  materials  exhibit  potential  advantages  over  noncomposite  materials 
because  they  can  display  unique  mechanical  prop>erties  and  characteristics  not 
possible  wiffi  a  single  constituent  Some  properties  potentiaUy  useful  in  marine 
applications  include  high  strength-  and  stiffiiess-to-weight  ratios,  relative  ease  of 
manufiicture  and,  in  comparison  to  ferrous  materials,  non-magnetic  property. 

Since  the  first  successful  application  of  fiber  reinforced  plastics  (FRP)  in 
the  1940's(^],  composite  materials  have  undergone  rapid  development.  During 
the  late  1960’s  and  early  1970's  advanced  composites  have  been  developed,  A 
new  group  of  plastic  laminates,  which  are  reinforced  by  fibers  (such  as  boron, 
graphic,  or  aramid-Kevlar  49)  with  higher  strength  and  stiffness  than  fiberglass, 
has  opened  up  a  new,  primary  structural  application  field  which  the  relatively 
low  modulus  fiberglass  composites  caimot  access. 

Yet,  because  of  the  intrinsic  complexities  of  mechanical  as  well  as  thermal 
interaction  among  constituent  materials,  the  theoretical  description  of  composite 
structures  has  advanced  much  slower  than  practical  applications.  Nevertheless  it 
was  appreciated  that  many  properties  of  FI^  have  proved  to  be  superior  to  those 
of  conventional  materials.  Consequently,  the  present  work  undertook  a 
preliminary  assessment  of  the  limitations  of  one  form  of  composite  structure  for 
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use  in  marine  propeller  blades.  This  assessment  intentionally  led  to  an 
examination  of  the  ability  of  present  structural  theory  to  support  critical  design 
decisions  involved  in  design  of  composite  propellers. 

The  problem  was  approached  by  performing  a  structural  evaluation  of  a 
typical  marine  propeller  blade  design,  which  was  constructed  of  a  hypothetical 
(thick-shell)  composite  material.  This  report  details  the  theoretical  analysis, 
hydrodynamic  loading,  and  finite-element  analysis  results.  The  resulting  stress 
and  deflection  behavior  is  discussed  and  compared  to  that  of  an  isotropic  metal 
blade  using  the  same  blade  form.  An  additional  theoretical  effort  was 
undertaken  to  develop  a  simplified  method  of  estimating  the  structural  behavior 
of  a  composite  blade,  which  did  not  require  extensive  effort  to  generate  finite- 
element-analysis  computer  input.  A  FORTRAN  code  is  presented  for  this 
method.  A  portion  of  the  material  in  this  report  was  previously  published  as 
Reference  4. 


STRUCTURAL  THEORY 

PROPERTIES  OF  A  LAMINA  (PLY) 

The  first  step  in  analyzing  a  composite  structure  is  the  determination  of 
elastic  properties  of  the  composite,  either  analytically  or  experimentally. 
However,  laboratory  testing  is  time  consuming  and  costly.  Thus  there  is  a  strong 
impetus  to  develop  an  analytical  tool  capable  of  predicting  composite  properties 
based  on  constituent  (fiber  and  matrix)  properties.  Though  in  the  past  many 
equations  were  derived  to  predict  the  mechanical  properties  of  composites,  most 
of  them  were  either  scattered  throughout  the  literature  or  not  of  simple  form. 
Usually  the  rule-of-mixtures  equations  were  used  to  smear  the  mechanical 
properties  of  fiber  and  the  matrix  into  effective  ply  (or  lamina)  properties. 
Chamis^^l  has  recently  developed  a  viable  formulation  for  the  elastic  modulus  of 
each  ply,  or  lamina,  by  the  use  of  composite  micromechanics.  His  unified  set  of 
equations  in  simple  form  is  foimd  to  useful. 

By  using  Chamis's  equations,  Nguyen  and  Critchfield^^J  have  presented 
an  analytical  capability  to  compute  the  complete  set  of  elastic  properties  for 
fibrous  composites.  The  following  is  a  brief  description  of  their  procedure : 

Elastic  Moduli  of  a  Lamina 

From  elementary  elasticity,  it  is  known  that  a  total  of  twenty-one 
constants  is  needed  to  describe  the  elastic  properties  of  an  anisotropic  material. 
For  orthotropic  materials  a  total  of  nine  independent  elastic  constants  will 
suffice.  If  the  material  of  each  lamina  is  assumed  to  be  transversely  isotropic, 
only  five  independent  elastic  constants  are  required.  These  elastic  constants 

(such  as  El ,  E2,  G12,  G23,  V12)  can  be  derived  by  using  constituent  material 
properties  following  Chamis's  method : 


2 
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El  —  Kf  Efi  +  Km  Em 

T-  Em  T- 

E2  = - = - =  E3 

l-VKf  (1-Em/Efi) 

Gi2  =  — J=Z — - =  Gi3 

1-V  Kf  (1-Gm/Gfi2) 

G23  = - . —  -  (1) 

1-V  Kf  (1-Gm/G£23) 


Vl2  =  KfVfi2+KmVm  =  Vi3 


2  , 
V23  = - 1 


where  En^  Ef2/  Qi2/  Gf23/  Vfi2  are  elastic  properties  for  the  individual  fibers  in 
various  directions;  Em/  Gm/  Vm  are  elastic  properties  for  the  matrix;  and  Kf,  i^n 
are  volume  fractions  of  fibers  and  matrix.  Ei  is  longitudinal  modulus;  E2  and  E3, 
transverse  moduli;  G12,  G13,  and  G23,  shear  moduli;  V12,  vi3^nd  V23  are  Poisson 
ratios.  Here,  Vfj  =  Poisson  ratio  for  transverse  strain  in  the  j-  direction  (ej)  over 
the  strain  (ei)  by  stressing  in  the  i-  direction. 

Roughly  speaking,  the  macroscopic  response  of  the  lamina  in  axial 
tension  is  assumed  to  be  predominantly  due  to  the  properties  of  fibers,  while  its 
behavior  in  both  transverse  tension  and  in  transverse  shear  will  be  mainly 
governed  by  the  properties  of  the  matrix. 

Elastic  Stiffness  of  a  Lamina 

Lamina  elastic  stiffnesses  based  on  the  predicted  lamina  moduli  can  be 
generated  from  the  lamina  stress-strain  relation 

(oi)  =  iQj]  {ej}  ;  i  ,j  =1,2,...,6  (2) 


where  Oj  =  stress  components 
Cj  =  strain  components 

Cij  =  stiffness  matrix  componenvS  or  lamina  elastic  stiffnesses 
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For  general  orthotropic  materials/  the  stifihess  matrix  can  be  constructed 
from  the  stress-strain  relation.  The  stiffness  matrix  components  for  a  layer  of 
orthotropic  material  are  given  in  Appendix  A. 

THIN-LAMINATE  (MULTIPLE-PLY)  COMPtDSITES 

A  laminate  is  an  assemblage  of  at  least  two  laminae  (thin  plies)  consisting 
of  parallel  fibers  stacked  (or  layed-up)  in  some  prescribed  manner  to  act  as  an 
integral  s^ctural  element.  Conventionally  the  layup  of  a  laminate  is  made  with 
the  ply  orientation  designated  by  the  angle  between  the  fiber  direction  and 
primary  load  direction  (see  Fig.  1). 


Fig  1.  An  off-axis  unicfirectional  lamina. 

For  general  laminates  where  the  layup  of  laminae  is  not  syirunetric  across 
the  thickness,  coupling  between  extension  and  bending  will  occur.  When  a 
laminate  is  subjected  to  external  loading,  the  resulting  forces  and  moments  bufit 
up  on  the  lammate  are  obtained  by  the  integration  of  the  stresses  in  each  layer  or 
lamina  through  the  laminate  thicimess.  The  integration  can  be  rearranged  to  take 
advantage  of  the  fact  that  the  stiffness  matrix  for  a  lamina  is  constant  within  the 
lamina,  and  with  the  substitution  of  the  lamiiui  stress-strain  relations,  one  can 
arrive  at  an  expression  relating  the  resulting  forces  (N)  and  moments  (M)  acting 
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on  a  laminate  to  the  middle  surface  strains  (£)  and  curvatures  (K)  of  the  laminate 
by  their  extensional  and  bending  stiffnessesJ^^^ 


N1  [  1(  e  )  .  .  2  6 

Ml  [  Bjj  Djj  JUr 


n 

where  Ay  «  (Qij)^  (r^-Zk-l ) 


8(1=1 1  (Qii)k  <4-^-1) 


Oil  =  Ui  +  U2  cos  26  +  U3  cos  46 


Qi2  -U4-U3  cos  46 
Qi6«^U2  sin  26  +  U3  sin  46 


Q22  =  Ui  -  U2  cos  26  +  U3  cos  46 
026  =  ^U2  sin26-U3Sin46 


Q66  =  U5  -  U3  cos  46 
TT  3Qii  +  3Q22  +  2Qi2  -»•  4Q66 


Q11-Q22 


,e  Q11  +  Q22-2Q12-4Q66 

U3 - g - 
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U4  = 


Oil  +  Q22  +  6Q12  -  4Q66 

8 


U5  = 


Oil  +  Q22-2Q12  +  4Q66 

8 


Qn  = 


El 

1-V12V21 


Qi2  = 


V12E2 

1-V12V21 


V21  El 
I-V12V21 


I-V12V21 


Q66  =  Gi2 

n  =  number  of  plies 

Z  =  vertical  (through-thethickness)  coordinate  of  each  ply  location 
(see  Fig.  2) 

It  should  be  noted  that  the  signs  in  the  expressions  of  Qis  and  Q26  in 
many  publications  (such  as  the  books  by  Jones^^i,  Vinson-Choul^i'  Nguyen  and 
Critchfield's  workl^l,  etc.)  are  incorrect.  They  should  be  corrected  as  shown 
above. 

THICK-LAMINATE  COMPOSITES 

It  has  been  often  shown  that  classical  bending  theory  without  the  effects 
of  transverse  shear  deformation  and  rotary  inertia  produces  errors  for  most  thin 
composite  laminatesi^l  Furthermore,  transverse  normal  stress  is  neglected  in  the 
classical  composite  laminate  theory.  Pagano  and  Pipes^^i  have  demonstrated 
that  high  interlaminar  normal  and  especially  shear  stresses  might  be  responsible 
for  onset  of  delamination  (splitting  at  the  interface  of  fiber  and  matrix)  near  the 
free  edge  of  the  laminate,  lire  effect  of  a  free  edge,  where  highly  localized  stress 
concentrations  often  occur  due  to  geometrical  as  well  as  materi^  discontinuities, 
is  a  critical  factor  in  the  theory  of  laminated  composites. 

In  view  of  these  deficiencies,  it  is  important  not  only  that  an  adequate 
analytical  tool  be  established  but  also  that  a  composite  be  so  designed  as  to 
provide  through-the-thickness  strength  to  eliminate  the  delamiruition.  These 
requirements  can  only  be  fulfiUed  by  a  threeKlimensional  analysis. 


DTRC/SHD-1373-01 


Fig.  2.  Geometry  of  an  n-layered  laminate. 

To  account  for  3-D  features  of  a  moderately  thick  laminate.  Sun  and 

Li  1^0]  presented  a  procedure  employing  a  long  wave  concept  to  derive  effective 
elastic  constants  for  thick  laminates  consisting  of  a  repeating  sublaminae.  They 
offered  two  sets  of  expressions: 

1 .  Full  expression — this  set  is  applied  to  general  cases  including  hybrid 
laminates; 

2.  Reduced  expression — ^this  set  is  the  simplified  version  of  the  full- 
expression  set  and  is  intended  for  a  single  composite  system. 

Theoretical  predictions  of  the  effective  elastic  constants  for  a  single 
composite  system  iMsed  on  the  reduced  expression  seem  to  check  very  weU  with 
some  other  niunerical  methods  and  even  with  some  limited  experimental 
results.’ 


Since  the  full  expression  formulation  is  of  lengthy  form  and  cannot  easily 
be  tmderstood,  an  alternate  formulation  based  on  the  same  assumptions  is 
proposed  as  follows : 


*  Private  communication  from  S.  Mayes,  DTRC  (1990) 
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For  a  sublaminate  containing  N  orthotropic  laminae  of  arbitrary 
thickness,  the  effective  elastic  constants  in  the  matrix  [  c  ]  can  be  obtained : 


N 

k=l 


Ak  =  C44<*'^  C45^^  -  (C45*y 

V|^  =  tj^  / 1  (tj^  is  the  thickness  of  the  kth  lamina,  and  t  is  the  total 
thickness  of  the  sublaminate) 

Cij^)  (i,  j  =  1, 2, 6)  :  stiffness  component  for  the  kth  ply ; 

For  a  single-composite  system,  the  above  formulation  reduces  to  that  of 
the  reduced  expression  given  by  Sun  and  Lil^^l. 
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WOVEN/NONWOVEN  FIBER  COMPOSITES 

3-D  textile  structural  composites^which  are  the  combination  of  a  resin 
system  with  a  textile  fiber,  yam  or  fabric  system  oriented  in  all  three  dimensions, 
have  recently  become  important  for  stmctural  applications!^^ I  Textile  stmctural 
composite  products  must  be  made  from  high-modulus  fibers  or  yams.  Yam  is  an 
assemblage  of  fiber  formed  into  a  continuous  strand  having  tex^e-like 
characteristics  distributed  periodically  in  the  nuiterial.  The  function  of  the  resin 
system  is  to  provide  rigidity  and  to  support  ar\d  protect  the  textile  reinforcement 
materials  (fibers,  yams,  or  fabrics)  in  a  prescribed  position  or  orientation  in  the 
composite. 

The  basic  fabric  sfructures  can  be  classified!^^!  as 

1 .  Wovens  —  made  by  interlacing  two  or  more  yam  systems  at  90* 
angles; 

2.  Krrits  —  made  by  interlooping  one  or  two  yams; 

3.  Braids  —  made  by  interwining  the  yam  system; 

4.  Nonwovens  —  made  by  mechanical  entanglement  of  fibers  or  by 
chemical  bonding  of  the  fibers  at  cross-over  points. 

The  comparison  of  basic  fabric  sfructures  is  described  in  Table  1  I 

Because  of  their  fully  integrated  nature,  such  as  fiber  continuity  and  3-D 
fiber  orientation  and  entanglement,  the  total  integrated,  advanced  fabric  systems 
are  thought  to  be  the  most  reliable  for  general  load-bearing  applications. 
Through-the-thickness  strength  and  interlaminar  strength  for  such  3-D  stmctural 
composites  have  been  greatly  increased  and  thus  the  risk  of  splitting  and 
delamination  is  minimized. 

The  engineering  design  and  analysis  of  3-D  fabric-reinforced  composites 
necessitates  an  analytical  capacity  that  can  relate  the  fabric  geometry  and 
processing  variables  to  the  mechanical  properties  of  the  composites.  However, 
very  little  is  known  regarding  this  relationship.  One  typical  procedure  based  on 
the  concept  of  unit  cell  geometry  was  established  by  some  investigators!^2-i6]  jq 
define  the  geometric  parameters  and  their  interrelationship.  For  example, 
Whyte!^^!  developed  a  Fabric  Geometry  Model  (FGM)  of  unit  cell  by  combining 
textile  engineering  methodology  and  a  modified  laminate  theory  to  demonstrate 
the  analysis  of  3-D  composites.  This  FGM  model,  specifically  for  3-D  braided 
metal  matrix  composite,  takes  into  consideration  the  material  properties,  fiber 
architecture,  and  processing  parameters.  By  using  the  fiber  volume  fraction 
semi-empirically  introduced  by  Ko!i21  together  with  the  properties  of  the 
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Table  1.  Comparison  of  basic  fabric  structures. 


Nonwoven 


Composition 

Formation 


Geometry 


yam 

interlace 


yam 

interloop 


QQQ 


Cell  model 


CW/  Imk 

Orientation 

Length 

Continuity 

Mobility 

Entanglement 


orthogonal 

short 

continuous 

limited 

flexible 


oblique 

long 

discontinuous 

tremendous 

mobile 


continuous 

limited 

flexible 


random 
very  short 
either 
very  slight 
rigid 


constituents,  the  stress-strain  relation  for  a  unidirectional  composite  can  be 
established  as: 

{o}  =  [c]{e}  ( 

The  stiffness  matrix  [cj  for  each  system  of  yams  can  be  calculated  by 
transforming  the  stiffness  matrix  [c]  of  a  unidirectional  composite  to  the 
appropriate  fiber  orientation  by  the  transformation  matrix  [Tj] : 

ICi]  =  [Ti][c]  [Til-1  < 

The  fabric-reinforced  composite  system  stiffness  is  determined  by 
superimposing  proportionaUy  the  stiffness  matrices  from  each  system  of  yams 
according  to  its  contributing  volume : 


[Csl  =  Ski[Ci] 
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where  [Csl  is  the  total  stiffness  matrix  and  kj  is  the  hractional  volume  of  the  ith 
system  of  yams. 

By  using  the  stiffness  matrix  [cj/  one  can  finally  determine  the  stress- 
strain  behavior  of  the  fabric-reinforced  composite. 

The  soundness  of  the  FGM  model  and  validity  of  the  results  for  3-D 
textile  stmctural  composites  need  further  theoretical  exploration  and 
experimental  verification. 

FAILURE  ANALYSIS 

There  are  many  different  ways  that  failure  can  occur  in  a  composite 
stmcture.  For  example,  the  possible  failure  modes  can  range  from  simple  loss  of 
structural  stiffness  due  to  yielding,  through  a  premature  warning  of  first-ply 
failure,  interlaminar  cracks,  or  delamination  by  fi-ee-edge  effect,  to  a  gross 
macroscopic  fracture.^^'^®'  While  fiiilure-analysis  work  has  been  substantially 
performed  on  composite  materials  in  the  past,  a  general  theoreticaUy  valid 
capability  does  not  exist.  Despite  this  lack,  normaUy  the  determination  of  failure 
modes  first  requires  obtaining  detailed  stress  (strain)  solutions.  To  achieve  this, 
one  often  has  to  rely  on  the  finite-element  method  as  no  other  method  can 
surpass  it  in  its  versatility.  However,  its  determination  of  stress  field  through  the 
assumption  of  a  simplified  displacement  field  in  the  classical  laiiunate  theory  is 
generally  not  reliable.  To  overcome  this  deficiency,  improved  and  efficient 
analytical  tools  are  required.  The  approach  by  Pagano^'^J  and  the  concept  by 
Sun  and  Maoi^^i,  all  referred  to  a  global-local  technique,  are  two  such  improved 
methods. 

Assuming  that  composite  stress  analysis  has  been  performed  and  that 
stresses  and  strains  are  available  everywhere,  one  then  can  simply  use  the 
maximum  stress  or  maximum  strain  theory  to  compare  the  predictive  values 
with  the  allowable  ones  in  each  particular  direction.  However,  most  composite 
materials  are  known  to  be  strongly  directionally  dependent.  With  no  interaction 
between  the  material  properties  in  different  directions  taken  into  account,  a 
failure  criterion  such  as  maximum  stress  or  strain  may  not  be  the  best  one  to  use 
when  a  composite  is  subjected  to  combined  stresses  or  strains.  The  most 
commonly  used  failure  criterion  seems  to  be  the  quadratic  failure  criterion  i2i-23l 
in  which  some  interaction  is  considered : 

Pxx  +  2  Fxy  Ox  Oy  +  Fyy  Oy2  +  Fjs  Og^  +  Fx  Ox  +  Fy  Oy  =  1  (11) 

where  Ox,  Oy  =  normal  stress  in  the  x,  y  directions  respectively 

Os  =  on-axis  shear  stress 
1 

^xx  =XX' 
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1 

Fyy  =yY' 

_i  J_ 

Fx  "X'X’ 

_i  J_ 

Fy  "Y'Y’ 

1 

Fss  =s2 

X(Y)  =  longitudinal  (transverse)  tensile  strength  in  the  x  (y) 
direction ; 

X'  (T)  =  longitudinal  (transverse)  compressive  strength  in  the 
X  (y )  direction ; 

S  =  longitudinal  shear  strength . 

Fxy  is  taken  to  be  (-1/2)^231  following  the  assumption  that  the  failure 
criterion  is  a  generalization  of  the  Von  Mises  criterion.  In  general,  equation  (11) 
should  be  m^ified  or  transformed  when  directions  of  applied  stresses  do  not 
coincide  with  the  material  axes. 

By  first  determiiung  the  maximum  length  of  the  strain  vectors  which 
satisfies  tiie  three-dimensional  ellipsoidal  strain  space  described  by  feilure 
criterion  equation  (11),  and  then  comparing  this  with  the  lengths  of  strain 
vectors  throughout  the  composite,  Daoust  and  Hoa^^Sl  are  able  to  determine  the 
safety  factor  in  composite  laminates  in  a  noniterative  fashion.  This  procedure 
and  some  other  theories  such  as  the  Tsai-Wu  theory ^^^1,  are  all  practical  methods 
for  predicting  composite  strength  since  test  data  for  combined  biaxial  and  shear 
strength  are  usually  difficult  and  expensive  to  obtain. 

SAMPLE  ANALYSIS  OF  A  PARTIAL-COMPOSITE  BLADE 
BLADE  GEOMETRY 

In  order  to  estimate  the  usefulness  of  available  theory  for  analyzing 
propeller  blade  structures,  the  stress  and  deflection  characteristics  of  a 
hypothetical  blade  were  predicted.  The  general  arrangement  of  the  blade 
configuration  is  shown  in  Fig.  3;  geometric  characteristics  are  listed  in  Table  2. 
The  blade  was  made  of  a  mixture  of  isotropic  and  composite  materials 
distributed  in  the  following  distinct  zones : 

The  tip  zone:  Type-4  nickel-aluminum  bronze  (NAB)  was  used  in  the 
blade  tip  area  from  xr  =  0.8  to  the  very  tip  end; 
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Fig.  3a  General  arrangement. 

Fig.  3.  Geometry  of  a  seven-bladed  thruster. 
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Table  2.  Geometric  characteristics  of  a  composite  thruster. 


(a)  Basic  Blade  Geomemc  paraineters 


Number  of  blade  (Z) 
Diameter  (D), 

Skew,  deg 
Hub-diameter  ratio 


7 

2 1.026  cm  (8.2779  in) 
30 

0.231 


r/R 

c/D 

9s  (dee.) 

it/d 

l/C 

f/fc 

.231 

.1718 

.5050 

-3.805 

-.0053 

.2138 

-.0243 

.3 

.1770 

.5561 

-7.402 

-.0114 

.1988 

-.0062 

.4 

.1830 

.7267 

-6.750 

-.0136 

.1772 

.014 

.5 

.1863 

.8857 

-2.458 

-.0061 

.1546 

.0285 

.6 

.1862 

.9345 

3.551 

.00922 

.1335 

.0333 

.7 

.1810 

.8966 

10.152 

.0253 

.1143 

.0317 

.8 

.1651 

.7849 

16.792 

.0366 

.0995 

.0242 

.9 

.1299 

.6367 

23.399 

.0414 

.0933 

.0069 

1.0 

.0695 

.0454 

29.998 

.0378 

.1075 

-.0488 

(b)  Blade  Surface  Form 


Blade 

Surface 

Offsets 

(in) 

r/R 

side 

.0 

.2 

Xc 

.4 

6 

.8 

1. 

.00 

.0960 

.1152 

.1150 

.0875 

.0134 

.231 

.00 

-.1403 

-.1815 

-.1813 

-.1317 

-.0134 

s 

.00 

.1075 

.1335 

.1333 

.0988 

.0116 

.3 

p 

.00 

-.1189 

-.1507 

-.1501 

-.1103 

-.0116 

s 

.00 

.1182 

.1517 

.1515 

.1100 

.0099 

p 

.00 

-.0905 

-.1101 

-.1099 

-.0823 

-.0099 

s 

.00 

.1205 

.1581 

.1579 

.1133 

.0088 

.5 

p 

.00 

-.0648 

-.0746 

-.0744 

-.0576 

-.0088 

s 

.00 

.1128 

.1496 

.1495 

.1067 

.0083 

.6 

p 

.00 

-.0471 

-.0511 

-.0510 

-0.411 

-.0083 

m 

s 

.00 

.0970 

.1291 

.1290 

.0922 

.OOD 

_ p _ 

.00 

-.0362 

-.0379 

-.0378 

-.0315 

-.0080 

o 

.00 

.0740 

.0981 

.0980 

.0707 

.0076 

.8 

IHhI 

.00 

-.0317 

-.0346 

-.0346 

-.0284 

-.0076 

■i 

.00 

.0438 

.0561 

.0561 

.0418 

.0073 

.9 

.00 

-.0342 

-.0418 

-.0418 

-.0322 

-.0073 

s 

.00 

.0061 

.0032 

0033 

.0051 

.0055 

1.0 

_E _ 

.00 

-.0420 

-.0571 

-.0572 

-.0411 

-.0055 

Note;  (i)  s  is  the  offset  from  nose-tail  line  to  suction  side; 

p  is  the  offset  from  nose-tail  line  to  pressure  side. 

(ii)  To  conven  (in)  to  (cm),  the  values  are  divided  by  0.3937. 


16 


DTRC/SHD- 1373-01 


The  composite  material  zone:  From  the  hub  to  xr  =0.8,  a  cross-section  of  this 
zone  reveals  a  sandwich  configuration  in  which  the  outer  thick-shell  skin  and  the 
internal  shear-web  support  were  made  of  fiber-reinforced  composites.  The  core 
is  filled  with  anisotropic  foam-type  material .  A  combined  total  of  approximately 
20  percent  of  the  chord  length  is  formed  of  low-modulus  materials  streamlined 
for  hydrodynamic  purposes  at  leading  and  trailing  edges. 

The  material  properties  used  in  the  different  zones  are  shown  in  Table  3. 
Finite  element  grid  lines  for  both  zones  at  the  isotropic-composite  juncture 
(xR  =0.8)  are  shown  in  Fig.  4. 

There  are  two  important  phases  in  dealing  with  the  3D  geometries  of  a 
complex  composite  blade.  First,  the  interior  boundary  between  skin  and  core  is 
mathematically  defined  and  3D  solid  elements  for  a  thick-shell  blade  are 
generated.  Such  meshing  of  solids  requires  the  accurate  repr^entation  of  the 
actual  blade,  fulfillment  of  mesh  continuity,  and  the  selection  of  well-shaped 
elements.  Second,  the  spatial  variation  or  orientation  of  material  properties 
should  be  determined  to  define  the  local  axis  system  for  material  definition. 

COMPOSITE  MATERIAL  CHARACTERISTICS 

The  analysis  of  a  composite  structure,  in  contrast  to  that  of  a  conventional 
single-material  structure,  is  such  that  elastic  material  constants  need  to  be 
determined  for  the  composite  system  embedded  in  the  structure.  The  analysis  of 
laminated  compxjsites  is  well  understood  and  can  be  readily  used.  However,  it 
has  been  shown  in  classical  composite  laminate  theory  that  the  onset  of 
delamination  often  occurs  near  the  free  edge  of  a  laminate,  where  highly 
localized  stresses  concentrate  due  to  geometric  and  material  discontinuities.  This 
free-edge  effect  can  be  minimized  or  practically  eliminated  by  using  2-D  or  3-D 
textile  structural  composites. 

A  particular  advantage  further  provided  by  3-D  braided  composites  is 
greatly  increased  through-the-thickness  or  interlaminar  strength,  reducing  the 
risk  of  delamination  Nonetheless,  3-D  composites  are  very  expensive  and 
the  determination  of  their  corresponding  elastic  constants  is  difficult  and  still 
state-of-the-art. 

In  the  present  study,  two-dimensional  braiding  fabric  was  chosen  as  the 
composite  material.  Braiding  fabric  can  be  readily  automated  and  manufiictured 
to  conform  to  complex  shapes  without  free-edge  effects.1241  For  predicting  elastic 
constants,  the  thick-laminate  technique^^^  was  adopted  to  account  for  3-D 
features,  along  with  the  assumption  that  composite  materials  are  layers  of 
orthotropic  materials.  It  was  further  assumed  that  the  skin  composite  byer  and 
shear-web,  consisting  of  brge  numbers  of  a  repeating  sublaminate,  could  be 
modelled  as  a  three-dimensional  homogeneous  anisotropic  solid.  Effective 
moduli  were  derived  by  properly  integrating  the  constituent  lamina  properties 
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Table  3.  Material  properties  of  a  composite  blade. 


Elastic  Moduli 

(msi)  and  Poisson 

Ratios  ( D ) 

Material 

Area 

Ell 

E22 

Gi2 

■U12 

013 

V23 

skin 

5.7 

1.4 

1.4  0.6  0.6 

0.36 

0.25 

0.25 

0.30 

composite 

web 

core 

same 

.04 

as 

.04 

above 

.04  .06  .06 

.016 

.25 

.25 

.25 

■I 

tip-2one 

edge 

E=18 

E  =  .01 

It  II 

p  P 

p  =  0.275  Ib/in^ 

P  =  0.037  Ib/in^ 

Note  :  (i)  Fiber  orientation  for  skin  (0/60/-60  degrees); 

fiber  orientation  for  shcar-web  (45/-45  degrees). 

Density  ( p  )  for  composite : 

0.073  Ib/in^  (  skin  and  web) .  o.020  Ib/in^  (core). 

(ii)  To  conven  (tnsi)  to  (N/m^),  the  values  are  multiplied  by  0.006894 

(iii)  To  convert  ( )  w  (  kg/m^ ).  the  values  are  multiplied 

by  0.00927. 
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through  the  thickness  of  the  sublaminate.  Three  different  kinds  of  composite 
materials,  as  shown  in  Fig.  5  and  Table  3,  were  used  in  the  study : 

1.  Monolithic  triaxial  E-glass  spar.  This  material  is  used  in  the  skin  layer 
of  the  composite  zone.  The  E-glass  fibers  were  oriented  at  three 
different  angles  (0, +60,-60  degrees)  to  the  element  axis  and  were 
woven  continuously  around  the  mandrel.  The  zero  angle  stands  for 
the  direction  of  unidirectional  fibers  running  parallel  to  the  braid  axis. 
Unlike  the  traditional  laminate  stack-up,  this  special  triaxial  shell  spar 
not  only  eliminates  the  free-edge  effect  but  also  is  suitable  for  taking 
up  the  major  bending  loads  produced  by  the  hydrodynamic  pressure. 

2.  ±45®  braided  E-glass:  This  material  was  used  in  the  shear-web  of  the 
composite  zone.  The  E-glass  fibers  were  oriented  at  ±45°  to  the 
element  axis.  This  ±45°  orientation  is  designed  to  resist  the  shear 
loads.  The  shear-web  also  increased  the  structural  stability  against 
buckling. 

3.  Polyurethane  core:  This  material  was  used  in  the  core  area  of  the 
sandwich  composite  zone.  The  core  material  can  not  only  resist  the 
shear  stress  produced  by  the  external  forces  but  also  has  the 
important  function  of  absorbing  the  energy  of  impacts. 

Assuming  orthotropic  elastic  symmetry  for  the  composite  materials,  only 
nine  independent  elastic  constants  for  each  of  the  different  composite  materials 
reqvured  for  the  3-D  finite  element  analysis  were  calculated  and  tabulated  in 
Table  4.  The  material  orientation  for  a  typical  composite  section  is  illustrated  in 
Fig.  6. 

BLADE  LOADING 

Two  kinds  of  blade  loads  were  considered  in  this  analysis; 
hydrodynamic  and  centrifugal.  Centrifugal  load  or  force  could  be  easily 
c^culated  since  it  is  a  function  only  of  blade  material  density,  blade  rotational 
speed,  and  blade  geometry.  Hydrodynamic  force  is  more  difficult  to  predict. 
There  are  several  available  prediction  codes.  The  PSF-2  numerical  code  was  used 
in  the  present  study  to  predict  the  pressure  distribution.  The  predicted  pressure 
coefficient  difference  generated  at  a  hypothetical  design  ship  speed  of  20  knots 
(10.29  m/s)  and  a  shaft  revolution  rate  of  n  =  4195  rpm,  applied  on  the  pressure 
side  of  the  blade,  is  plotted  in  Fig.  7  and  tabulated  in  Table  5. 
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URETHANE  ACRYLATE  RESIN 


MONOLITHIC  TRIAXIAL  E-GLASS  SHELL 


±45°  BRAIDED  E-GLASS  SHEAR  WEB 


POLYURETHANE  FOAM  CORE 


URETHANE  ACRYLATE  RESIN 

Fig.  S.  Material  arrangement  in  composite  zone. 

Table  4.  Elastic  constants  of  composite  blade  for  finite-element  analysis. 


Elasdc  Constants  (msi) 

for  Finite  -Element  Analysis 

Area 

Dim 

Dn22 

D2222 

Dn33 

D2233 

D3333 

D1212 

Di313 

D2323 

skin 

3.26 

1.04 

3.26 

0.515 

0.515 

1.58 

1.11 

0.45 

0.45 

web 

2.75 

1.55 

2.75 

0.515 

0.515 

1.58 

1.62 

0.45 

0.45 

core 

.048 

.016 

.048 

.016 

.016 

.048 

.016 

.06 

.06 

Note  :  (i)  For  the  definitions  of  stiffness  element  Dim  etc,  one  is  referred  to  the 
ABAQUS  manual. 

(ii)  To  conven  (msi)  to  (N/m^),  ihc  values  arc  multiplied  by  0.006894. 
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3  •  through-the-thlcknes8  direction 


Fig.  6.  Orientation  for  materiai  definitions. 
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- i — i —  (  psf2  -100  )  /dcp  :  xr  :  0.2650 

...Z.-i-..  <  psf2  -100  )  /dcp  :  xr  r  0.3250 

- s — 5—  (  psf2  -100  )  /dcp  :  xr  r  0.3756 


Fig.  7.  Pressure  differerK:e  (^p)  versus  chordwise  station  (xc) 
for  different  stations  (psf-2  code). 
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Table  5.  Pressure  difference  on  composite  blade. 


Pressure 

Difference 

ACd 

r/R 

.025 

.25 

.35 

Xc 

.45 

.55 

.65 

.75 

.85 

.975 

.265 

-2.675 

-.451 

-.287 

-.108 

-.062 

-.099 

-.104 

-.109 

-.194 

.325 

-3.825 

-.402 

-.118 

.199 

.319 

.308 

.338 

.358 

.068 

.375 

-3.934 

-.176 

.167 

.542 

.700 

.705 

.736 

.725 

.217 

.425 

-3.584 

.148 

.525 

.924 

1.104 

1.117 

1.131 

1.064 

.331 

.475 

-3.066 

.506 

.903 

1.307 

1.495 

1.511 

1.499 

1.369 

.431 

.550 

-2.472 

.985 

1.404 

1.812 

2.006 

2.014 

1.961 

1.750 

.584 

.650 

-2.091 

1.442 

1.896 

2.322 

2.522 

2.512 

2.411 

2.118 

.674 

.750 

-2.062 

1.561 

2.037 

2.476 

2,677 

2.660 

2.553 

2.262 

.734 

.825 

-1.873 

1.261 

1.680 

2.068 

2.252 

2.255 

2.197 

1.999 

.727 

.875 

-1.369 

.768 

1.058 

1.330 

1.463 

1.478 

1.463 

1.361 

.546 

.925 

-.349 

-.035 

.014 

.062 

.089 

.102 

.117 

.128 

.095 

.975 

.684 

-.607 

-.778 

-.935 

-1.013 

-1.021 

-1.011 

-.947 

-.365 

FINITE-ELEMENT  STRESS  ANALYSIS 
ABACUS  Computer  Code 

From  the  geometric  input  and  the  material  orientation  at  the  centroid  of 
each  solid  element  in  the  composite  zone,  the  finite  element  code  ABAQUS^^Sl 
was  used  to  perform  the  stress  analysis.  The  ABAQUS  code  is  a  general  purpose 
finite-element  computer  code  for  production  use  in  a  wide  range  of  applications. 
It  has  an  extensive  library  of  nonlinear  features.  Input  is  organized  in  key  words 
and  sets.  The  ABAQUS  code  has  been  successfully  applied  in  many  engineering 
problems.  For  example,  the  tip  motion  of  a  geometrically  nonlinear  isotropic 
cantilever  beam  with  a  circular  pipe  cross-section  was  obtained  by  using 
ABAQUS  (see  example  1.2.1  in  Ref.  24)  and  compared  excellently  with  the  exact 
solution  for  the  inextensive  beaml^l.  An  anisotropic  composite  plate  was 
analyzed  by  using  ABAQUS(see  example  4.4.1  in  Ref.  25}  and  yielded  deflection 
and  moment  values  which  compared  closely  with  the  analytical  solution.  1271 

Two  different  kinds  of  elements  were  used  in  the  blade:  20-node  and  15- 
node  curved  solid  elements.  They  were  adopted  to  represent  the  curve-  and 
skew-geometry  of  the  entire  blade.  There  were  a  total  of  430  solid  elements. 
Except  for  six  15-node  triangular  prisms  used  in  the  high  curvature  leading  edge 
of  the  composite  zone,  the  remaining  are  20-node  elements.  The  boimdary 
condition  required  that  nodal  points  on  the  intersection  of  blade  and  hub  be  held 
fixed.  Care  was  taken  along  the  xr  =  0.8  section,  the  dividing  juncture  between 
isotropic  and  composite  zones,  that  MFC  (multi-point  constraints)  were  properly 
employed  to  ensure  structural  connectivity. 
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Failure  Criterion 


No  failure  criterion  tailored  especially  for  composite  materials  was 
employed  for  this  analysis.  Instead,  to  simplify  the  analysis,  stresses  were 
considered  to  exceed  allowable  levels  when  the  elastic  limit  of  the  material  was 
exceeded.  For  this  purpose  the  effective  stresses  were  obtained  from  the  effective 
modulus  model  of  Sun  and  (In  a  more  general  treatment,  failure  criteria 
would  be  applied  to  individual  through-the-skin  stress  levels  recovered  from 
effective  stresses  computed  by  the  global-local  technique^^l) 

Numerical  Results 

Ntunerical  results  were  obtained  for  the  blade  with  a  partial  composite- 
shell  structure,  and,  for  comparison,  a  second  blade  of  solid  NAB,  which,  of 
course,  is  isotropic.  The  same  geometry  and  loading  were  used  in  both  cases. 

Calculated  elastic  deflections  for  both  the  composite  and  isotropic  blades 
are  shown  in  Figs.  8  and  9.  As  can  be  seen,  the  deflection  pattern  of  the 
composite  blade  is  similar  to  that  of  the  isotropic  blade.  However,  the  nuiximum 
deflection  at  the  tip  of  the  composite  blade,  listed  in  Table  6,  is  an  order  of 
magnitude  larger  than  that  of  the  isotropic  blade.  This  relatively  large  deflection 
has  the  potential  of  causing  a  significant  change  in  geometric  shape. 

Calculated  stresses  in  the  radial  direction  on  the  suction  and  pressure 
sides  for  the  composite  blade  (excluding  the  isotropic  tip  zone)  are  plotted  in 
Figs.  10  and  11.  The  principal  stresses  for  the  isotropic  blade  are  plotted  in 
Bgs.  12  and  13.  These  principal  stress  values  are  predominantly  in  the  radial 
direction,  and  therefore  correspond  closely  to  radial  stress  values.  Shear  stresses 
on  the  pressure  side  for  composite  and  isotropic  blades  are  plotted  in  Hgs.  14 
and  15,  respectively.  Through-the-thickness  tensile  stress  for  the  composite 
blade  is  plotted  in  Fig  16. 

By  comparing  Fig.  10  with  Fig.  12,  it  is  seen  that  the  stress  contours  on  the 
suction  side  for  both  composite  and  isotropic  blades  exhibit  a  similar  trend,  in 
that  the  maximum  stress  is  located  near  the  mid-chord  of  the  blade-hub 
intersection,  with  smaller  stress  magnitudes  toward  the  tip  and  leading  edge. 
Similarly,  Figs.  11  and  13  show  that  stress  contours  on  the  pressure  side  also 
exhibit  a  similar  pattern  for  both  composite  and  isotropic  blades,  with  the 
maximum  stress  occurring  near  the  trailing  edge  (0.85  chord-length  measured 
from  leading  edge)  aroxmd  the  radial  station  xr=0.35.  Maximum  shear  stresses 
(Figs.  14  and  15)  also  occur  close  to  the  trailing  edge  near  the  root 

A  more  detailed  examination  was  made  of  the  through-the-thickness 
tensile  stress  distribution.  The  maximum  through-the-thickness  tensile  stress  is 
seen  from  Fig.  16  to  occur  on  the  inner  surface  of  the  thick-shell  skin  layer  at  the 
trailing  edge  near  the  root.  The  stress  distribution  in  the  thickness  direction  was 
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Table  6.  Maximum  stress  and  deflection  values. 


NAB  Blade  •• 

Composite  Blade** 

Response 

FEA 

Simple 

Be^ 

FEA 

Simplified 

Method 

Suss 

-5060 

-5900 

-8350 

-7060 

Stress* 

Slips 

5916 

7130 

8050 

7490 

(psi) 

S12 

2280 

2250 

3067 

1528 

S33 

— 

— 

2600 

— 

Deflection 

S 

0.0073 

— 

0.072 

0.052 

(in)  I  I _ I _ 

*  Sii  :  Bending  stress  in  the  radial  direction  (for  simple  beam 


this  is  referred  to  principal  stress; 

$12  :  Shear  stress  on  the  blade  cross-section  in  the  skin  layer 
parallel  to  the  blade  surface; 

S33  :  Through-the-thickness  tensile  stress  in  the  skin  layer 
thickness  direction; 

(+)  Tensile  stress;  (-)  Compressive  stress; 

ss :  suction  side;  ps :  pressure  side;  5 :  elastic  deflection  at  blade 

tip. 

•*  The  working  stress  levels  of  NAB  are  taken  to  be 
approximately*^) 

511  (tension)  =  9000  psi;  Sn  (compression)  =  -90(X)  psi; 

512  =  4500  psi.  For  composite  triaxial  E-glass  material,  ultimate 
strength  stress  values  were  estimated  by  the  quadratic  failure 
criterion(211  to  be; 

511  (tension)  =  8100  psi;  Sn  (compression)  s  -356(X)  psi; 

512  =  6300  psi.  For  S33,  matrix  ultimate  strength  is  expected  to  be 
about  3000  psi. 

Note:  To  convert  (psi)  to  (N/m^),  the  values  are  multiplied  by 
6894;  to  convert  (in.)  to  (cm),  the  values  are  divided  by  0.3937. 
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Fig.  10.  Compressive  radial  stress  contour  on  suction  side  (composite  blade). 
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Fig.  11 .  Tensile  radial  stress  contour  on  pressure  side  (composite  blade). 
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Fig.  12.  Principal  (compressive  stress  contour  on  suction  side  (solid  blade). 


Fig.  13.  Principal  (tensile)  stress  contour  on  pressure  side  (solid  blade). 
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Fig.  16a.  Blade  planform. 


Fig.  16.  Contours  of  through-the-thickness  tensile  stress 
on  the  pressure  side  of  composite  blade. 
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such  that  the  maximum  tensile  stress  generated  on  the  inner  shell  surface 
decreased  to  a  smaller  stress  level  toward  the  outer  skin  layer. 

The  maximum  stresses  and  deflections  for  both  NAB  and  composite 
blades  are  summarized  in  Table  6,  along  with  NAB  blade  stress  values  calculated 
by  a  simple  beam  theory  customarily  used  for  propeller  design  work.  For  NAB, 
the  maximum  stresses  predicted  by  FEA  did  not  exceed  66  percent  of  the 
allowable  working  levels  as  indicated  in  the  table  notes.  Beam  theory 
consistently  gave  higher  values,  up  to  79  percent  of  the  allowable  levels.  Both 
results  imply  that  the  NAB  blade  design  would  be  adequate  for  operation  at 
higher  than  the  hypothetical  design  speed. 

The  composite  stresses  were  compared  against  estimated  ultimate 
strength  levels.  The  in-plane  predicted  stresses  reached  over  99  percent  of  the 
ultimate  stress  levels,  implying  that  the  blade  had  reached  its  failure  point  at  the 
design  speed  used.  Similarly,  interlaminar  stresses  reached  87  percent  of  an 
approximate  ultimate  stress  for  r  matrix.  These  results  indicate  that  the 
composite  blade  would  be  limited  to  operation  at  speeds  well  below  the  design 
value  used  in  the  analysis. 

In  evaluating  these  results,  it  must  be  noted  that  the  weight  savings  for 
the  composite  blade  are  more  than  50  percent  when  compared  to  the 
conventional  blade  made  of  metal  alone. 


SIMPLIHED  STRESS/DEFLECTION  CALCULATION 
FOR  A  COMPOSITE  BLADE 

It  is  usually  quite  involved  and  time  consuming  to  predict  the  detailed 
stresses  and  deflections  of  a  composite  blade  using  a  3D  finite  element  method. 
To  expedite  the  calculation  process  during  the  preliminary  design  phase,  it  is 
often  desirable  and  useful  for  a  designer  to  have  a  simplified  method  to  quickly 
estimate  the  conventional  bending  and  shearing  stresses  and  deflections.  In 
what  follows,  a  simple  scheme  for  predicting  stress  and  deflection  is  briefly 
described. 

DETERMINATION  OF  FORCES  AND  MOMENTS 

In  a  right-handed  Cartesian  reference  frame,  the  unit  base  vectors  are 
( i ,  j  ,  k  )  for  the  x-,  y-,  and  z-axes.  The  x-axis  is  positive  downstream;  y-axis  is 

positive  starboard;  z-axis  is  positive  upward.  The  hydrodynamic  force  (Fj, )  at 
each  radial  cross  section  is 

Eh  =-Fh“*1>pi  +  Ph“"1’p'e 
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where  £9  =  -  cos  0  j  -  sin  0  k 
0  =  tan-i  (-y/z) 

The  centrifugal  force  (F^ )  at  each  section  is 
F  =F  e 

i,c  c  ilr 

# 


where  e  =  -  sin  0  j  +  cos  0  k 

The  corresponding  hydrodynamic  moment  (Mh)  and  centrifugal  moment 
(Me)  are 

Mh  =  dxF, 

M  =  dxF 


where  d  is  the  displacement  vector  between  each  section. 

Assuming  that  the  root  section  is  fixed,  the  blade  acts  like  a  cantilever 
beam.  The  moments  at  each  section  are  determined  by  integrating  the 
corresponding  forces  from  tip  (free)  end  down  to  the  section  of  interest. 

CALCULATION  OF  BLADE  STRESS 

The  bending  stresses  are  approximately  calculated  as  follows: 

g|,(ps)l 

<^^(88)  j  \  \ 

where  =  In-plane  moment  about  the  nose-tail  line 
tk  =  Maximum  blade  cross  section  thickness 

Ik  =  Moment  of  inertia  about  the  nose-tail  line 
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aic  =  Cross  section  area  occupied  by  composite  skin  layer 

Sk  =  Shell  skin  layer  thickness  in  the  composite  blade  section 

k  =  Radial  section  station  indec 

ps  =  pressure  side 

ss  =  suction  side 

The  shear  stress  is  calculated  as 

_ _ 

where  Qk=  Torque  in  the  radial  direction 
Ak  =  Total  blade  cross  section  area 

Assuming  that  the  root  section  is  fbced,  the  blade  acts  like  a  cantilever 
beam.  The  moments  at  each  section  are  determined  by  integrating  the 
corresponding  forces  from  tip  (free)  and  down  to  the  section  of  interest. 

CALCULATION  OF  TIP  DEFLECTIONS 

In  this  simplified  deflection  calculation,  only  the  in-plane 
moments  and  the  corresponding  bending  rigidities  about  the  nose-tail  line  of 
each  section  were  considered.  The  bending  rigidity  of  each  section  was 
calculated  as  the  product  of  the  effective  elastic  modulus  (following  [10])  in  the 
bending  axis  and  the  moment  of  inertia  of  the  cross  section  about  the  nose-tail 
line,  for  which  each  section  was  represented  as  an  elliptical  shell  having  its  major 
axis  coincident  with  the  noise-tail  line.  This  approximation  introduces 
inaccuracy  because  the  noise-tail  line  might  not  pass  through  the  neutral  axis, 
and  the  section  is  not,  in  general,  elliptical.  Improved  accuracy  can  be  obtained 
by  externally  determining  the  moment  of  inertia  and  inserting  those  values  into 
the  computation.. 

Once  the  distribution  of  in-plane  moments  (M(x))  and  the  bending 
rigidities  (EI(x))  are  known,  the  maximum  tip  deflection  is  computed  as 


M(x)  (R-x) 


EI(x) 


dx . 
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To  evaluate  this  simplihed  method,  the  stress  and  deflection  of  the 
previously  analyzed  composite  blade  were  calculated.  The  program  input  and 
resulting  output  are  given  in  Appendix  B,  along  with  a  FORTRAN  computer 
code  for  these  calculations. 

RESULTS  OF  SIMPLIHED  CALCULATION 

Stress  values  were  calculated  at  nine  stations  from  the  root  to  the  tip,  for 
the  mid-chord  location.  Maximum  stresses  occurred  at  the  root  section.  The 
maximum  stress  values,  which  are  given  in  Table  6,  ranged  from  7  to  50  percent 
below  the  finite-element  values. 

The  deflection  calculation  also  showed  a  lower  value  with  the  simplified 
approach.  As  shown  in  Table  6,  the  calculated  tip  deflection  was  28  percent 
below  the  finite  element  value. 


DISCUSSION 


ANALYTICAL  APPROACHES 

The  flnite-element  approach  proved  suitable  for  relatively  comprehensive 
analysis  of  the  performance  of  the  composite  blade.  The  treatment  level 
appeared  sufficient  for  use  in  assessing  the  adequacy  of  a  particular  blade  design 
for  at  least  preliminary  design  purposes.  It  is  true,  however,  that  not  all  aspects 
of  the  composite  structure  were  represented  in  the  present  analysis. 
Simplifications  were  made  because  of  the  extraordinary  complexity  of  composite 
structures,  and  because  further  work  would  have  required  a  disproportionate 
investment  of  resources. 

Composite  blade  design  would  benefit  from  an  even  simpler  structural 
design  approach.  However,  as  illustrated  by  the  present  simplified  method, 
further  simplification  can  result  in  the  loss  of  significant  features  (such  as 
interlaminar  stress)  of  the  structure.  A  further  apparent  deficiency  of  the  present 
approach  was  its  imconservative  stress  predictions,  unlike  the  beam  analysis  for 
a  solid,  isotropic  blade,  which  gave  conservative  results.  These  results  show  that 
further  exploration  is  needed  of  approaches  for  simplified  stress  and  deflection 
predictions. 

Shell  structure  performance  typically  is  governed  by  bending  stresses  in 
the  outer  layer,  and  by  shear  stiffness  of  webs  and  core.  Composite  shell  analysis 
must  pay  attention  to  through-the-thickness  stress  as  well,  because  of  the 
relatively  weaker  strength  of  the  matrix  which  may  predomiiwte  in  the  through- 
the-thiclmess  direction.  This  finding  suggests  that  simplified  approaches  may 
not  be  able  to  be  based  on  two-dimensional  stress  models. 
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Even  if  stress  distributions  could  be  determined  reliably,  application  of 
meaningful  failure  criteria  remains  as  a  problem.  A  primary  deficiency  of  the 
present  analyses,  and  perhaps  of  the  held  in  general,  is  absence  of  detailed 
failure  criteria.  In  composites,  failure  processes,  such  as  fiber  rupture  and 
buckling,  fiber /matrix  debonding,  matrix  yielding/cracking,  and  delamination, 
can  occur  in  a  very  complex  way.  Furthermore,  many  factors,  such  as  the 
residual  stresses  due  to  fabrication,  nonlinear  and  inelastic  effects,  and  interface 
interactions  between  fibers  and  matrix,  do  not  appear  to  be  properly  accounted 
for  in  available  theory. 

As  a  result  of  these  deficiencies,  current  approaches  to  composite  material 
strength  and  failure  are  semi-empirical  in  nature.  These  approaches  must  be 
developed  to  a  sufficient  degree  to  support  propeller  blade  design  if  such  blade 
structures  are  to  be  accepted  for  vehicular  use.  Simplified  methods,  while  highly 
desirable  for  preliminary  design  estimates,  may  be  problematical.  The  large 
deviations  found  between  the  simplified  and  finite-element  methods  suggests 
that  a  large  error  allowance  would  be  required  if  the  present  simplified  method 
were  to  be  used  for  preliminary  design. 

The  author  would  like  to  emphasize  the  need  for  development  of  fetigue 
limits  for  evaluating  the  feasibility  of  blades  such  as  the  present  one.  Such  limits 
would  presumably  include  techniques  for  experimentally  certifying  candidate 
structures,  such  as  are  used  for  current  metal  blades. 

COMPOSITE  BLADE  PERFORMANCE 

The  present  blade  example  was  intended  to  demonstrate  the  similarities 
and  differences  between  conventional  NAB  and  composite  thick-shell  blade 
structures.  Substitution  of  a  lower-modulus  FRP  material,  with  the  resulting 
much-larger  stress  and  deflection  values,  did  not  produce  an  equivalently 
performing  propeller  blade.  Rather,  it  provides  a  frame  of  reference  to  assist 
designers  who  are  just  beginning  to  explore  applications  of  the  new  material  to 
become  fiimiliar  with  what  may  be  radically  different  features  and  parameter 
values.  In  this  vein,  both  structural  and  hydrodynamic  designers  may  need  to 
modify  their  designs  and  analysis  methods  to  exploit  potential  benefits  of 
composite  materials. 

One  such  new  behavioral  characteristic  of  composite  blades  may  be  larger 
deflections.  Acceptability  of  a  large  tip  deflection,  in  the  present  case  about  1.7 
percent  of  the  overall  propeller  radius  and  an  order  of  magnitude  greater  than 
that  of  the  NAB  blade,  must  be  evaluated  in  terms  of  powering  performance.  No 
study  of  the  hydrodynamic  implications  of  the  deflection  was  performed. 

A  second  new  characteristic  is  stress  levels  for  which  no  established 
interpretation  exists.  Stress  values  will  require  interpretation  in  terms  of 
established  working  levels.  Even  without  established  levels,  of  course,  the 
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present  stress  levels  were  clearly  unacceptable  at  the  speed  used.  At  a  minimum, 
the  present  composite  blade  would  be  limited  to  speeds  well  below  those  at 
which  the  NAB  blade  would  reach  its  working  stress. 

A  third  new  characteristic  is  interlaminar  stresses.  Such  stresses  may 
cause  failure  at  levels  which  would  be  acceptable  for  intralaminar  stresses.  The 
present  example  suggests  that  a  critical  region  for  such  stresses  may  be  the  small- 
radius  aft  edge  of  the  blade  strength  member. 

Despite  the  above  differences,  there  was  considerable  similarity  between 
the  qualitative  stress  patterns  in  the  plane  of  the  composite  shell  and  parallel  to 
the  solid  blade  surface.  In  that  respect,  the  composite  blade  may  exhibit 
relatively  "conventional"  behavior. 

Lastly,  a  complete  feasibility  study  of  a  composite  blade  structure  should 
include  variations  in  material  properties  such  as  fiber  stiffness  and  strength  as 
well  as  blade  section  geometry,  which  were  not  included  in  the  present  study. 

CONCLUSIONS 

1 .  Three-dimensional  finite-element  stress  analysis  appears  to  be 
adequate  for  preliminary  design  of  composite  blades. 

2.  The  present  composite-adapted  b^m  theory  approach  for  rapidly 
estimating  stresses  and  deflections  in  thick-sh^  composite  blades  is 
not  adequate  for  use  in  preliminary  blade  design. 

3.  Establishment  of  acceptable  working  stress  levels  for  composite 
propeller  blade  structures  is  required  for  successful  blade  design. 

4.  Use  of  composite  blade  structures  may  be  possible  in  more 
applications  if  larger  blade  deflections  than  for  conventional  metal 
blades  can  be  accepted. 

RECOMMENDATIONS 

1.  A  satisfactory  method  for  rapidly  estimating  blade  stresses  in 
preliminary  design  should !:«  developed. 

2.  Working  stress  levels  should  be  determined  for  composite-material 
propeller  blades. 

3.  Effects  of  large  blade  deflections  on  blade  hydrodynamic  performance 
should  be  determined. 
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APPENDIX  A. 

STIFFNESS  MATRIX  FOR  GENERAL  ORTHOTROPIC  MATERIALS 

For  general  orthotropic  materials,  the  stiffness  matrix  is  given  by  the 
following  stress-strain  relation : 


Cll  Ci2  Ci3  0  0  0 

Ci2  C22  C23  0  0  0 

C13C23C33  0  0  0 

0  0  0  C44  0  0 

0  0  0  0  C55  0 

_  0  0  0  0  0  C66 


where 


<^1,  <^2,  ^3  =  normal  stresses  in  the  1, 2, 3  directions ; 

Ti2>  '^23,  'C31  =  shearing  stresses  in  the  1-2, 2-3,  and  3-1  planes ; 

61  >  £2,  £3  =  normal  strains  in  the  1, 2, 3  directions 

Yl2>  Y23,  Y31  =  shearing  strains  in  the  1-2, 2-3,  and  3-1  planes. 

Note  that  there  is  no  interaction  between  normal  stress  and  shearing 
strain  or  shearing  stress  and  normal  strain  and  that  there  are  only  nine 
independent  constants  in  this  sti^ess  matrix. 

The  stiffness  matrix  components  for  a  layer  of  orthotropic  material  are; 


Cii  = 


1  -  V23  V32 
E2  E3  A 


Ci2  = 


V21  +  V31  V23  V12+V13  V32 


E2  E3  A 


El  E3  A 


Cl3  = 


V31  +  V2I  V32  Vi3  +  V12  V23 


E2  E3  A 


El  E2  A 
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C22  = 


1  -V13V3I 
El  E3  A 


C23  = 


V32  +  V12V3I 
El  E3  A 


V23+V21  Vl3 
El  E2  A 


C33  = 


1  -  V12  V21 
El  E2  A 


C44  =  Gi2 

C55  =  G23 

C66  =  G31 

El,  E2,  E3  =  Young's  moduli  in  1, 2,  and  3  directions,  respectively 


Vjj  =  Poisson's  ratio  for  transverse  strain  in  the  j-direction  (Ej)  over  the 
strain  (Ei)  by  stressing  in  the  i-direction 

where 


1 .  V12V21-V23V32-V31  Vi3-2Vi2  V32V13 
El  E2  E3 
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APPENDIX  B. 

PROGRAM  BSINP  —  A  SIMPLIFIED  METHOD  OF  STRESS  AND 
DEFLECTION  CALCULATIONS  FOR  COMPOSITE  BLADES 

INPUT  FORMAT 
The  input  to  BSIMP  is  read  as  the  following 
read  32,title 

read  kd,  nxr,  dim,  rpm,  sthk,  eskin,  eprop 

read  *,  (yr(i),  i=nxr,l,-l) 

read  *,  (ph(i),  i=nxr,l,-l) 

read  *,  (rake(i),  i=nxr,l,-l) 

read  (th(i),  i=nxr,l,-l) 

read  (c(i),  i=nxr,l,-l) 

read  (t(i),  i=nxr,l,-l) 

read  (fdp(i),  i=nxr,l,-l) 

read  *,  (wt(i),  i=nxr,l,-l) 

read  *,  (saKi),  i=nxr,l,-l) 

read  *,  (seix(i),  i=nxr,l,-l) 

read  *,  (peix(i),  i=nxr,l,-l) 

The  variables  are  defined  below. 

Card  1:  Title 

This  is  a  72-character  title  for  the  problem  identification. 

Card  2:  kd,  nxr,  dim,  rpm,  sthk,  eskin,  eprop 

kd  is  an  integer  variable  identifying  the  the  type  of  imits  used  in 
the  input. 


DTRC/SHD-1373-01 


43 


kd  =0  indicates  that  the  unit  will  be  of  dimensional  form; 
kd  >0  indicates  that  the  unit  will  be  of  non-dimensional  form; 
nxr  Number  of  input  radii; 
dim:  Thruster  diameter; 

rpm:  Thruster  shafter  rotational  speed  in  revolution  per 
minute; 

sthk:  Shell  skin-layer  thickness  in  the  composite  blade  section; 

eskin:  Effective  elastic  modulus  of  skin-layer  in  the  bending 
axis; 

eprop:  Elastic  modulus  of  solid  blade  with  isotropic  material. 
Card  3:  (yr(i),  i=nxr,l/-l) 

This  card  contains  mcr  values  of  nondimensional  radius 
card  4:  (ph(i),  i=nxr,l,-l) 

If  kd  =0,  this  card  contains  nxr  values  of  blade  pitch  angle  in 
degrees; 

If  kd  >0,  this  card  contains  nxr  values  of  pitch-to-diameter  ratio. 

Card  5:  (rake(i),  i=nxr,l,-l) 

If  kd  =0,  this  card  contains  nxr  values  of  total  rake; 

If  kd  >0,  this  card  contains  nxr  values  of  total-  rake-to-diameter 
ratio. 

Card  6:  (th(i),  i=nxr,l,-l) 

This  card  contains  nxr  values  of  skew  angle  in  degrees. 

Card  7:  (c(i),  i=nxr,l,-l) 

If  kd  =0,  this  card  contains  nxr  values  of  chord-length; 

If  kd  >0,  this  card  contains  nxr  values  of  chord-to-diameter 
ratio. 
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Card  8:  (t(i),  i=nxr,l,-l) 


If  kd  =0,  this  card  contains  nxr  values  of  blade  section 
maximum  thickness; 

If  kd  >0,  this  card  contains  nxr  values  of  thickness-to-chord 
ratio. 

Card  9:  (fdp(i),  i=nxr,l,-l) 

This  card  contains  nxr  values  of  hydrodynamic  force  per  unit 
length. 

Card  10:  (wt(i),  i=nxr,l,-l) 

This  card  contains  nxr  values  of  centrifugal  force  |3er  unit 

length. 

Card  11:  (sar(i),  i=nxr,l,-l) 

This  card  contains  nxr  values  of  cross-sectional  area  of 
composite  blade. 

Card  12:  (seix(i),  i=nxr,l,-l) 

This  card  contains  nxr  values  of  effective  bending  rigidities  (=EI) 
of  composite  blade  section. 

Card  13:  (peix(i),  i=nxr,l,-l) 

This  card  contains  nxr  values  of  bending  rigidities  (=Ei)  of 
isotropic,  solid  blade  section. 


VARIABLE  DEHNITIONS 
The  major  outputs  are  the  following  variables: 
xr  Nondimensional  radius. 

fip:  In-plane  force  in  the  nose-tail  line  (positive  toward  the  trailing 
edge). 

fop:  Out-of-plane  force  (positive  toward  suction  side), 
fri:  Force  in  the  radial  direction. 
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mip:  In-plane  moment  about  the  nose-tail  line  (positive  toward  the 
trailing  edge). 

mop:  (Dut-of-plane  moment  (positive  toward  suction  side), 
mz:  Torque  in  the  radial  direction, 
smps:  Tensile  stress  at  midchord  of  blade  pressure  side, 
smss:  Compressive  stress  at  midchord  of  blade  suction  side, 
tau:  Shearing  stress  at  midchord  close  to  the  blade  surfece. 

EXAMPLE  INPUT  FOR  COMPOSITE  BLADE 

0  9  8.2779  4195.0  0.04  0.285  1.8A 

0231  0.3  0.4  05  0.6  0.7  0.8  0.9  1.0 

34.833  30.542  30.040  29.417  26.371  22.181 17.344  12.691  8223 

-0.0441  -0.09465  -0.11279  -0.05006  0.07630  020929  030307  034258  031316 

-3.805  -7.402  -6.750  -2.458  3.551  10.152 16.792  23399  30.0 

1.42225  1.46519 1.51486 1.54202  154150 1.49830  136689  1.07561  0.57531 

03038  02914  0.2682  0.2384  0.2061  0.1714  0.1358  0.1002  0.0621 

0.0000  -3.4010  1.7283  11.0722  20.0709  24.4383  205959  5.0384 

0.0000 

0.755E-02  0.769E-02  0.781E-02  0.780E-02  0.765E-02  0.729E-02 

0.653E-02  0502E-02  0.256E-02 

0.103  0.105  0.107  0.107  0.105  0.999E-01  0.894E-01 

0.688E-01  0350E-01 

0.348E-03  0324E-03  0.276E-03  0211E-03  0.147E-03  0.902E-04 

0.449E-04  0.150E-04  0.195E-05 

0353E-02  0320E-02  0.259E-02  0.185E-02  0.119E-02  0.665E-03 

0304E-03  0.956E-04  0.120E-04 
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EXAMPLE  OUTPUT  FOR  A  COMPOSITE  BLADE 


kd=  0  ( Dimensional  input ) 


kd,ns,dim(in),rpm,sthk(in),eskin(10**7  psi),eprop(10**7  psi)  =  0  9  8.278 

4195.000  0.040  0.285  1.800 


XT  =  0.2310 

0.9000  1.000 

0.3000 

0.4000 

0.5000 

0.6000 

0.7000 

0.8000 

ph  (deg)  =  34.83 

12.69  8.223 

30.54 

30.04 

29.42 

26.37 

22.18 

17.34 

rake  (in)  =-0.4410E-01  -0.9465E-01  -0.1128  -0.5006E-01 

0.3031  0.3426  0.3132 

0.7630E-01  0.2093 

th  (deg)  =  -3.805 
23.40  30.00 

-7.402 

-6.750 

-2.458 

3.551 

10.15 

16.79 

c  (in)  =  1.422 

1.076  0.5753 

1.465 

1.515 

1.542 

1.542 

1.498 

1.367 

t  (in)  =  0.3038 

0.1002  0.6210E-01 

0.2914 

0.2682 

0.2384 

0.2061 

0.1714 

0.1358 

fdp  (Ib/in)  =  0.0000 

-3.401 

1.728 

11.07 

20.07 

24.44 

20.60 

5.038  0.0000 


wt  ab/in)  =  0.7550E-02  0.7690E-02  0.7810E-02  0.7800E-02  0.7650E-02 
0.7290E-02  0.6530E-02  0.5020E-02  0.2560E-02 


sar(sq.in)  =  0.1030  0.1050  0.1070  0.1070  0.1050 

0.8940E-01  0.6880E-01  0.3500E-01 


0.9990E-01 


seix(10**7psi)=  0.3480E-03  0.3240E-03  0.2760E-03  0.2110E-03  0.1470E-03 
0.9020E-04  0.4490E-04  0.1500E-04  0.1950E-05 


peix(10*^psi)=  0.3530E-02  0.3200E-02  0.2590E-02  0.1850E-02  0.1190E-02 
0.6650E-03  0.3040E-03  0.9560E-04  0.1200E-04 
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i 

xr 

fip 

fop 

fri 

mip 

mop 

mz 

1 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

2 

0.900 

0.123 

1.070 

3.021 

0.064 

0.559 

-0.117 

3 

0.800 

0.491 

6.547 

7.088 

0.629 

2.448 

-1.375 

4 

0.700 

0.553 

16.267 

11.017 

4.214 

5.654 

-5.927 

5 

0.600 

0.302 

26.111 

14.091 

11.362 

9.829 

-14.332 

6 

0.500 

0.092 

33.317 

16.199 

21.415 

13.786 

-24.584 

7 

0.400 

0.820 

36.610 

17.839 

33.818 

15.456 

-32.316 

8 

0.300 

0.683 

36.380 

19.943 

48.121 

15.414 

-32.976 

9 

0.231 

-3.097 

35.085 

22.279 

58.487 

18.958 

-26.542 

- •••••  Deflection  at  tips  ***** - 

deflection  (in.)  solid/hollow  =  0.5608E-02  0.5159E-01 


• 

1 

XT 

smps 

smss 

tau 

1 

0.231 

7492.080 

-7059.476 

1158.489 

2 

0.300 

6357.235 

-5977.375 

1464.808 

3 

0.400 

4849.579 

-4516.136 

1528.174 

4 

0.500 

3599.255 

-3296.462 

1312.994 

5 

0.600 

2404.279 

-2135.874 

914.579 

6 

0.700 

1251.428 

-1030.875 

492.272 

7 

0.800 

350.336 

-191.765 

172.170 

8 

0.900 

-17.326 

105.136 

29.749 

9 

1.000 

0.000 

0.000 

0.000 

Fortran  STOP 

48 


DTRC/SHD- 1373-01 


LISTING  OF  COMPUTER  CODE  BSIMP 


program  bsimp 
implicit  real*8  (a-h,o-z) 
parameter  (nx=15) 
character*72,title 

dimension  x(nx),y(nx),2(nx),th(nx)/rake{nx),rx(nx),ph(nx)^m(nx), 
1  c(nx),wt(nx),fdp(nx),cf(nx),cfy(nx)/:fe(nx),fip(nx),fop(nx), 
2£r(mt)/bmip(m(),bmop(nx),tmz(nx),bmih(nx),bmoh(nx),tmh(nx), 
3  bmic(nx),bmoc(nx),tmc(nx)pcr(nx),bml  (nx),bm2(nx)^ebc(nx), 
4peix(nx)^i(nx),pei(nx),smei(nx),pmei(nx),yr(nx)^bc(nx),pix(nx) 
dimension  tmq(nx)/ir(nx)/ii(nx),t(nx),tk(nx)^r(nx)^(nx), 

1  sigps(nx)^igss(nxXtau(nx) 

c —  kd  =0  (dim  lb-in  unit),  kd  >0  (nondim  unit:p/d4t/d,c/d,t/c,f/c) 
pi=3.141592654 
rd=pi/180.0 
jl=l 

read  32,title 

read  *,kd,nxr/Jin\,rpm,sthk,eskin,eprop 

read  *,(yr(i),i=:nxr,l,-l) 

read  *,(ph(i),i=nxr,l,-l) 

read  *,(rake(i),i=nxr,l,-l) 

read  *,(th(i),i=nxr,l,-l) 

read  *,(c(i),i=nxr,l,-l) 

read  •,(t(i),i=nxr,l,-l) 

read  •,(fdp(i),i=nxr,l,-l) 

read  *,(wt(i),i=nxr,l,-l) 

read  *,(sar(i),i=nxr,l,-l) 

read  •,(seix(i),i=nxr,l,-l) 

read  •,(peix(i),i=nxr,l,-l) 

print  M,' - output  of  ',title 

if  (kd  .eq.  0)  print  31, 'kd  =  ’,kd,'  ( Dimensional  input )’ 
if  (kd  .ne.  0)  print  31,'kd  =  ’,kd,'  ( Nondimensional  input )' 
if  (kd  .eq.  0)  go  to  10 
do  11  i=l,nxr 

ph(i)=datan(ph(i)/yr(i)/pi)/rd 

rake(i)=rake(i)'^dim 

c(i)=c(i)*dim 

t(i)=t(i)'^c(i) 

11  continue 
10  continue 

print  24,'kd,ns/lun(in),rpm,sthk(in),eskin(10*'*7  psi),eprop(10*'V 

1  psi)  =’,kd,nxr,dim,rpm,sthk,eskin,eprop 

print  23,'xr  =',(yr(i),i=nxr,l,-l) 

print  23,'ph  (deg)  =',(ph(i),i=nxr,l,-l) 

print  23,'Fake  (in)  =',(rake(i),i=nxr,l,-l) 

print  23,'th  (deg)  =',(th(i),i=nxr,l,-l) 
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print  23/c  (in)  =’,(c(i),i=nxr,l,-l) 

print  23,’t  (in)  =',(t(i),i=nxr,l  ,-l ) 

print  23,'fdp  (Ib/in)  =',(fdp(i),i=nxr,l,-l) 
print  23,’wt  Ob/in)  =',(wt(i),i=nxr,l,-l) 
print  23/sar  (sq.in)  =',(sar(i),i=nxr,l/-l) 
print  23/seix(10**7  psi)  =',(seix(i),i=nxr,l/-l) 
print  23/peix(10**7  psi)  =',(peix(i),i=nxr,l,-l) 
rad=0.5*din\ 
iad23:rad*rad 
print  26 

romga=rpm*pi/30.0 
fop(l)=0.0 
fip(l)=0.0 
bmip(l)=0.0 
bmop(l)=0.0 
tm2(l)=0.0 
fxi=0.0 
fyi=0.0 
^=0.0 
bmx=0.0 
bmy=0.0 
bm2=0.0 
bmip(l)=0. 
bmop(l)=0. 
tniz(l)=0.0 
fxhisO.O 
fyhiaO.O 
fehi=0.0 
fyci=0.0 
fed=0.0 
bhx=0.0 
bhysO.O 
bh2=0.0 
bcx=0.0 
bcy=0.0 
bcz=0.0 
rx(l)=yr(l)*rad 
do  1  i=2,nxr 
rx(i)=yr(i)*rad 
ths=th(i)*rd 
thsml=th(i-l)*rd 
x(i)=rake(i) 
x(i-l)=rake(i-l) 
sth=dsin(ths) 
cth=dcos(ths) 
sthml  =dsin(thsml ) 
cthml  =dcos(th5m1 ) 


50 


DTRC/SHD- 1373-01 


y(i)=-rx(i)*sth 

z(i)=  rx(i)*cth 

y  (i-1  )=-rx(i-l  )*sthml 

z(i-l)=  rx(i-l)*cthinl 

dx=x(i-l)-x(i) 

dy=yO-l)-y(i) 

dz=z(i-l)-z(i) 

dr=rx(i-l)-rx(i) 

adr=dabs(dr) 

sm(i)=wt(i)/32.2 

cf(i)=sm(i)*(rx(i)/  12.0)*romga*rom^ 
cfy(i)=-cf(i)*sth 
cfz(i)=  cf(i)*cth 
sm(i-l)=wt(i-l)/32.2 

cf (i-1  )=sin(i-l  )*(rx(i-l )  / 1 2.0)*roin^*romga 

cfy  (i-1  )=-cf (i-1  )*sthml 

c£z(i-l)=  cf(i-l)*cthinl 

phi=ph(i)*rd 

sphi=dsin(phi) 

q>hi=dcos(phi) 

phiml  =ph(i-l  )*rd 

sphiml  =dsin(phiml ) 

q)himl =dcos(phiml ) 

C]fKt=q)hi*cth 

q>st*qphi*sth 

spct=sphi*cth 

spst=sphi*sth 

qxtml =q)hiinl  ’cthml 

q>stinl  =q)himl  •sthml 

spctml =sphiml  *cthml 

spstml  ssphiml  •sthml 

fhxu=-0.5*fdp(i-1  )*qjhiml*adr 

flixb=-0.5*fdp(i)*q)hi*adr 

fhyu=-0.5*fdp(i-l  )*spctml  ’adr 

fhyb=-0.5*fdp(i)*spct*adr 

fhzu=-0.5*fdp(i-l  )*spstml  •adr 

£hzb=-0.5*fdp(i)*spst*adr 

fcyu=0.5*cfy(i-l  )*adr 

fcyb=0.5*cfy(i)*adr 

fczu=0.5*cfe(i-l  )*adr 

fczb=0.5*cfe(i)*adr 

fhxi=fhxu+fhxb 

fhyi=fhyu+fhyb 

fhk=fh^+fhzb 

ftyi=fcyu+fcyb 

fczisfczu+fczb 

bdml=fxi 
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fyiml=fyi 

feimlsfzi 

fxhiml=£xhi 

fyhiml=fyhi 

fehiinl=fehi 

fyciml=fyd 

£zciinl=fizci 

fxhi=fxhi+fhxi 

fyhi=fyhi+fhyi 

fzhi=fehi-i-fhzi 

fyci=fyci+fcyi 

^=fzd+fczi 

fxi=fcci+fhxi 

fyi=fyi+fhyi+fcyi 

W=£zi+fhzi+fczi 

fip(i)=fxi*sphi-fyi*q)ct-fei*q3st 

fop(i)=-fxi*q)hi-fyi*spct-fei*spst 

fr(i)=-fyi*sth+fei*cth 

bmhx=(0.75*fhzu+0.25*fhzb)*dy-(0.75*fhyu+0.25*fhyb)*dz 

bmhy=(0.75*fhxu+0.25*fhxb)*dz-(0.75*fhzu+0.25*fhzb)*dx 

bmhz=(0.75*fhyu+0.25*fhyb)*dx-(0.75*fhxu+0.25*fhxb)*dy 

bincx=(0.75*fczu+0.25*fczb)*dy-(0.75*fqru+0.25*fcyb)*dz 

bmcy=-(0.75*fczu+0.25*fczb)*dx 

bmcz=  (0.75*fcyu+0.25*fcyb)*dx 

if  (i  .i\€.  2 )  go  to  3 

bmx=bmx+bmhx+bmcx 

bmy=bmy+bmhy+bincy 

bmz=binz+bmhz+bmcz 

bhx=bhx+bmhx 

bhy=bhy+bmhy 

bhz=bhz-(-bmhz 

bcx=bcx+bmcx 

bcy=bcy+bmcy 

bcz=bcz+bmcz 

go  to  4 

3  bmx=binx+bmhx+bmcx+fziml*dy-fyiml*dz 
bmy=bmy+binhy+bmcy+fximl  Mz-feiml  Mx 
bmz=binz+bmhz+bmcz+fyiml  *dx-fximl  *dy 
bhx=bhx4-binhx+fzhiin1  •dy-fyhiml  *dz 
bhy=bhy+bmhy+fxhiml  •dz-fzhiml  *dx 
bhz=bhz+bmhz+fyhiml  •dx-fxhiml  *dy 
bcx=bcx+bmcDc+feciinl  *dy-fyciinl  •dz 
bcy=bcy+bmcy-fzciinl  ’dx 
bcz=bcz+bincz+fyciml  *dx 

4  bmip(i)=binx*sphi-bmy*cpct-bmz*cpst 
bmop(i)=-bmx*cphi-bmy*spct-bmz*spst 
tinz(i)=-bmy*sth+bmz*cth 
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bmih(i)=bhx*sphi-bhy*cpct-bhz*qjst 

bmoh(i)=-bhx*q)hi-bhy*spct-bhz*spst 

tmh(i)=-bhy*sth+bhz*cth 

bmic(i)=bcx*sphi-bcy*cpct-bcz*cpst 

bmoc(i)=-bcx*q)hi-bcy*spct-bcz*spst 

tmc(i)=-bcy*sth+bcz*cth 

1  continue 
c  print  25 

c  do  5  i=24ucr 

c  print  214,bnuh(i),bmoh(iXtmh(i),bmic(i),bn\oc(i),tinc(i) 
c  5  continue 
print  20 
do  2  i=l,nxr 

print  21,i,yr(i),fip(i),fop(i),fr(i),bmip(i),bmop(i),tmz(i) 

2  continue 
binip(l)=0.0 
do  7  i=l,nxr 
ii=nxr-i+l 
xr(i)=yr(u) 
bnnl(i)=bmip(ii) 
tmq(i)=tmz(ii) 
cfr(i)=fr(ii) 
ch(i)=c(ii) 
tk(i)=t(ii) 
sa(i)=sar(ii) 
six(i)=seix(ii)/eskin 
pix(i)=peix(ii)  /eprop 
sei(i)=seix(ii)*10.0*'7 
pei(i)=peix(ii)*l  0.0**7 

7  continue 
print  26 

c  print  23/  xr  =',(xr(i),i=  1  ,nxr) 

c  print  23/  bm  =',(bnil  (i),i=l  ,nxr) 
c  print  23/  seix  =',(seix(i),i=nxr,l,-l ) 

c  print  23/  peix  =',(peix(i),i=nxr,l,-l) 

c  print  23/  sei  =',(sei(i),i=l,nxr) 

c  print  23/  pei  =',(pei(i),i=l,nxr) 

do  6  i=l,nxr 
xrl=1.0-xr(i) 

smei(i)=bm1  (i)*xrl  /sei(i) 
pmei(i)=bml(i)*xrl  /pei(i) 
c  print  1 9,i,xrl  ,bml  (i)4ei(i)^mei(i),pmei(i) 

6  continue 

c  print  23/  smei  =',(smei(i),i=l4'xr) 

c  print  23,'  pmei  =',(pmei(i),i=l,nxr) 

print  27 

defc=l  .0'^rad2*simpun(xr(l  ),smei(l),nxr) 
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defp=l  .0*rad2*siinpun(xr(l  ),pmei(l  ),nxr) 

print  23/  deflection  (in.)  solid/hollow  =  ’/iefp4efc 

print  26 

do  8  i=l,nxr 

ty=0.5*tk(i) 

ax=0.5*(ch(i)-sthk) 

by=ty-0.5*sthk 

az=pi*ax*by 

yps=ty 

yss=-ty 

sigps(i)=bm1  (i)*yps/six(i)+cMi)/sa(i) 
sigss(i)sbm1  (i)*yss/six(i)+cfr(i)  /sa(i) 
tau(i)=-0.5*tmq(i)/(az*sthk) 

8  continue 
print  29 

do  9  i=l,nxr 

print  214,xr(i),sigps(i),sigss(i),tau(i) 

9  continue 

19  format  (i3,12(g9.2,lx)) 

20 formate/'  i  xr  fip  fop  fri 
1  mip  mop  mz ') 

21  format  (i3,7(fl2.3,2x)) 

22  format  (2i3,12(fl0.3,lx)) 

23  format  (a,10(gl2.4,lx)) 

24  format  (a^4^,10(f8.3,lx)) 

25 formate/'  i  bmih  bmoh  tmh  bmi 

Ic  bmoc  tme  ') 

26  format  (' - ') 

27  format  (/' - *****  Deflection  at  tipjs  ***** - '/) 

28  format  (/) 

29 formate  i  xr  smps  smss  tau') 

30  format  (/a^/) 

31  format  (a,i3^,a) 

32  format  (a) 
stop 

end 

FUNCTION  SIMPUNe  X,  Y,  N ) 

IMPLICIT  REAL’S  (A-H,0-Z) 

C 

C  FORTRAN  IV  FUNCTION  FOR  SIMPSON’S  RULE  INTEGRATION 
C  ARBITRARY  NO.  AND  LENGTH  INTERVALS 
C 

DIMENSION  X(15),  Y(15) 

C 

IF(N.LE.1)THEN 

C 

C  N  .LE.  1  -  •’’-►ERROR’’*- 
C 
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non  non  non  non  nn 


RETURN  ZERO  ANSWER 


SIMPUN  =  0.0 
ELSEIF(N.EQ.2)THEN 

N  =  2  -  USE  TRAPAZOIDAL  RULE 

SIMPUN  =  (Y(l)  +  Y(2))*(X(2)  -  X(l))/2.0 
ELSE 

N  .GE.  3  -  USE  SIMPSON'S  RULE 
IF  ( MOCK  N,  2 )  .EQ.  0 )  THEN 

N  IS  EVEN,  TAKE  A  SPECIAL  STEP  TO  START  THE  INTEGRATION 

DX21  =  X(2)-X(1) 

DX31  =  X(3)-X(1) 

DX32  =  X(3)  -  X{2) 

S  =  DX21*(  Y(l)*(3.0  -  DX21/DX31)  + 

1  Y(2)*(3.0  +  DX21  /DX32)  - 

2  Y(3)*DX21'“2/(DX31*DX32)  )/6.0 
L*3 

ELSE 

N  IS  ODD,  NOTHING  SPEQAL  REQUIRED 

S  =  0.0 
L  =  2 
END  IF 

DO  1000  K  =  L,  N-1, 2 
ABSDX  =  dABS(  X(K)  -  X(l) ) 

IF  ( dABS(  X(K-l)  -  X(l) )  .GE.  ABSDX  .OR. 

1  ABSDX  .GE.  dABS(  X(K+1 )  -  X(l) ) )  GO  TO  1100 

DXKPKM  =  X(K+1)  -  X(K-1) 

DXKKM  =  X(K)-X(K-1) 

DXKPK  =  X(K+1)-X(K) 

S  =  S  +  DXKPKM*(  Y(K-1P(3.0  -  DXKPKM/DXKKM)  + 

1  YCKPd  .0  +  DXKPKM/DXKKM  +  DXKKM/DXKPK)  + 

2  Y(K+ir(2.0  -  DXKKM/DXKPK)  )/6.0 
1000  CONTINUE 

SIMPUN  =  S 
END  IF 
GO  TO  99999 
1100  CONTINUE 
print  100,  K,X(K) 

SIMPUN  »  0.0 
99999  CONTINUE 
RETURN 

100  FORMAT(23HONON  MONOTONE  X  SIMPUN  ,I4,1PE12.4) 

END 
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